Adiabatic invariance with first integrals of motion 
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The construction of a microthermodynamic formalism for isolated systems based on the concept 
of adiabatic invariance is an old but seldom appreciated effort in the literature, dating back at least 
to P. Hertz [Ann. Phys. (Leipzig) 33, 225 (1910)]. An apparently independent extension of such 
formalism for systems bearing additional first integrals of motion was recently proposed by Hans H. 
Rugh [Phys. Rev. E 64, 055101 (2001)], establishing the concept of adiabatic invariance even in 
such singular cases. After some remarks in connection with the formalism pioneered by Hertz, it will 
be suggested that such an extension can incidentally explain the success of a dynamical method for 
computing the entropy of classical interacting fluids, at least in some potential applications where 
the presence of additional first integrals cannot be ignored. 
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Consider an isolated collection of particles with / de- 
grees of freedom evolving in time according to the usual 
laws of classical particle mechanics (e.g., a gas of N = 
//3 particles inside a closed cylinder). The dynamics of 
this system can be uniquely described by a trajectory in 
the 2/-dimensional phase space of the canonically con- 
jugate variables q = (<Zi, (//) and p = (pi, ...,pf). Let 
-ff (q, p; A) be the Hamiltonian of the system depending 
(for simplicity) on a single external parameter A (e.g., the 
volume set by the position of the piston in the cylinder). 
Assume further that no other first integrals exist (see be- 
low, however). Then, for constant A, the trajectory of 
the system will be confined to and will (almost) fill out 
the whole surface of constant energy H(q, p; A) = Ex, 
provided ergodicity is assumed. Consider now a pro- 
cess in which A = X(t) varies very slowly in compari- 
son to a typically small "observation time" r, which in 
turn is a sufficiently large quantity with regard to micro- 
scopic processes such that time averages of phase space 
functions taken over it approximate the corresponding 
microcanonical averages at constant A (this is an ideal- 
ized instance of an otherwise precisely observed process, 
e.g. a reversible compression/expansion of the gas in the 
cylinder above). Such processes are usually referred to 
as adiabatic (though an additional quasistatic qualifier 
would be certainly very welcome), and it was shown by 
P. Hertz in his papers titled "On the Mechanical Foun- 
dations of Thermodynamics" ^] that the (r-averaged) 
phase space space volume Q(E, A) enclosed by the sur- 
faces of constant energy H(q, p; A) = E\ is a conserved 
quantity during processes of this kind (justifying then 
the label adiabatic invariant) , much like the macroscopic 
entropy of Clausius under the same conditions. Ther- 
modynamics could thus be constructed from mechanical 
arguments, leading directly to a microscopic form of en- 
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tropy, Sq = lnf2, which should be contrasted to the usual 
form inspired by Boltzmann's ideas involving the phase 
space area u> — dfl/dE, S w = In u, the difference between 
these two being particularly relevant for finite systems 
H f§ j§. These results of Hertz were greatly appreciated 
by Einstein (see fl in particular) and are often taken as 
the starting point of statistical mechanics in the German 
literature J?], H. However, these simple and elegant ob- 
servations seem to have escaped the attention of modern 
H [l(], |ll| and even classic |l^] treatments (some rare 
exceptions are Refs. Jl3|, fl4|). The underlying derivation 
is simple and accessible and thus will not be reproduced 
here (this so-called ergodic adiabatic invariance problem, 
however, is far from being exhausted by Hertz's origi- 
nal papers, being constantly addressed at different levels 
in the literature, see e.g. (15| [U| ]l§|). It suffices 
to mention here that, in connecting Rugh's derivation 
with Hertz's results, an infinitesimal change of the pa- 
rameter A under the presence of parameter-independent 
first integrals Fj(q, p) = Ii maintains straightforwardly 
the adiabatic conservation of the phase space volume 
£l(E, I, A), with / = {Ii}, and thence of the entropy Sq, 
this observation following immediately from the fact that 
(dfl/dli) dli = under a change dX (compare, for exam- 
ple, Eq. (I 236) of Ref. § or Eq. (34.4) in Ref. §). An- 
other important result presented in Ref. Jigf] was the pos- 
sibility of computing the "bulk" temperature by means 
of a microcanonical (and therefore temporal) average at 
constant E, I and A, namely T n = { Y • VH(Y)\E, I, A), 
with Y parallel to the surfaces Fj(q, p) = Ii and satisfy- 
ing V • Y = 1. I note in passing that a similar form of 
this generalized equipartition theorem was obtained by 
Miinster Q for the case of arbitrary cyclic coordinates 
present in the Hamiltonian (which when cast in terms of 
appropriate variables, e.g. center-of-mass for the conser- 
vation of linear and angular momentum or action-angle 
for more general integrals, translates essentially into the 
above requirement of parallelism). 

Finally, I would like to bring attention to the fact that 
the extension provided by Rugh of Hertz's original ideas 
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can incidentally explain the success of the so-called "adi- 
abatic switching method" Q for computing entropies of 
classical interacting fluids (which is one of the very few 
works to make explicit use of Hertz's results), at least 
for the applications discussed below. This method takes 
full advantage of the adiabatic invariance of Sa by ini- 
tially considering a "reference" system with Hamiltonian 
Hq whose entropy is known explicitly (e.g., an ideal gas) 
and slowly turning on the interactions such that the final 
Hamiltonian equals the desired one, Hi. If this process is 
sufficiently slow, it can be considered as adiabatic (in the 
sense defined above), the phase space volume Q(E, A) 
(and thus the entropy Sq) being a conserved quantity 
throughout the whole switching process p4|] . Therefore, 
a thermodynamic quantity which is not a phase space 
function and hence not immediately obtained by means 
of the usual microcanonical/temporal average [plf can 
be in principle easily computed. However, this method 
is usually applied in the "molecular-dynamics ensemble" 
in which the total linear momentum P is exactly con- 
served due to the use of periodic boundary conditions 
(see, e.g., pi], sec. 2.10] and references therein). Such ad- 
ditional first integrals of motion should be carefully con- 
sidered when using the Hertz invariant, since the man- 



ifold accessible to the system clearly does not coincide 
with the whole constant-energy surface. This fact seems 
to have passed unnoticed by the authors of Ref. [ p0[ , 
although the above observations by Rugh really give a 
strong theoretical support to their method. Indeed, it is 
easy to see that all three components of P are parameter- 
independent in this case: the switching on of the inter- 
nal interactions (which depend on the relative position 
of the particles only) clearly does not break the transla- 
tional invariance of the Lagrangian. This same idea can 
be equally well applied to a more correct formulation of 
this ensemble that only recently was realized |22|, |23| , in 
which three additional first integrals (the components of 
the center of mass R) are considered. It is worth empha- 
sizing that this simple result was only possible because 
of the parameter-independence of the integrals of motion. 
The opposite and more general case, however, is still an 
open issue, as pointed out by Rugh. 
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